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PART I 

Introduubion. 

In o.i’der to determine the feasibility of the space shuttle radar it 
was necessaiy to implement a computer simulation program at the Uni- 
versify of California, Berkeley. This program simulates the spectrum 
which results when a radar signalls transmitted into the ionosphere for 
a finite time and received for an equal finite interval. The Spectrum 
derived from this signal is statistical in nature because the signal is 
scattered from the ionosphere, which is statistical in nature. Thus, any 
estimator which uses the backscattered spectrum has a statistical char- 
acter. That is, many estimates of aiy properly of the ionosphere can 
be made. Their average value will approach the average properly of 
the ionosphere which is being measured. I-Iowever, due to the statistical 
nature of the spectrum itself, the estimators will vary about this average. 
The square root of the variance about this average is called the standard 
deviation, an estimate of the error which exists in any particular radar 
measurement. In order to de ter mine the feasibility of the space shuttle 
radar, the magnitude of these errors for measurements of physical 
interest must be understood. 

A complete reading of tliis report makes it possible to implement 
the computer simulation described herein. If such an implementation 
is not desired, reading of Part II will give a good understanding of the 
physical principles involved in implementation of the program and a 
reasonable overview of the mathematical tools used . 
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PART n 

The ore iical Considerations and Computer Operation 

The computer program simulates a spectrum having statistical pro- 
perties which are estimated in frequency space. Before explaining 
exactly how the computer performs this simulation, it is necessary to 
discuss the physical ionosphere that is being simulated. This ionosphere 
is time stationaiy. This means that at any particular point in the iono- 
sphere being simulated, a long time average will yield a well defined 
number for any physical properly being measured. It also means that 
an ensemble average of many ionospheres will yield the same results 
as a long time average of ary particular ionosphere being studied. The 

signals back to the radar. This does not mean that there are no statistical 
fluctuations within the ionosphere, but, rather, that the point pro per ties 
measured anywhere in the ionosphere possess the same long time average . 
Thus, although at any particular time there may be density fluctuations 
in the ionosphere, on the average the density at any point in the ionosphere 
is identical with tlie density of any other point from which scattered 
signal is being returned to the radar. The real ionosphere is not uniform. 

It is often true that over the time period during which data from the radar 
are being averaged or over the spatial region that one pulse of the radar 
covers, the ionosphere varies in its properties. This type of non-uniformity 
is not considered in this stu(fy, and the uniformity assumption should be 
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considered in extrapolating the data from the simulation to real situations 
in the ionosphere, 

In addition to the uniformity of the scattering medium, the computer 
simulation is idealized in another respect. It is assumed throughout 
this study that the radar itself is ideal. For example, no asymmetries 
of the received filters are simulated. Nor is it assumed the receiver 
gain varies with time due to saturation hy the transmitted signal, although 
both of these problems exist in real radars. 

Consider a time stationary ionosphere in which densily fluctuations 
can exist, and compare it with an ionosphere in which no density fluc- 
tuations exist. Assume that a signal with wavelength, X, long by com- 
parison to the Debye length and frequency, co, high by comparison to the 
plasma frequency is transmitted into a fluctuationless ionosphere. There 
will be no return backscattered Signal into the radar receiver. This can 
be shown as follows. Consider two planes in the ionosphere perpendicular 
to the vector of the incident radar wave. For convenience call the two 
planes PI and P2 and allow the first of these to be at a distance d from 
the radar and the secohd at a distance d + A /4. During the trip from 
radar to scattering plane and back, the scattered signal from the first 
plane travels a distance, 2d, while the scattered signal from plane P2 
travels a distance 2d +x/2. The signals from these two planes are 
out of phase with one another and will cancel at the radar receiver. Now 
consider any plane, P3, in the medium from which scattered signal is 
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expected. It is possible to pick another plane, P4, which is related 
to P3 in the same way that PI and P2 were related. In this way it can 
be seen that no scattered signal can return to the radar if the ionosphere 
is fluctuationless. 

Now consider the time stationary ionosphere in which fluctuations 
exist. Consider a plane PI and another plane P2 as before, but which 
have larger densities than the SurroundiJig medium. These planes are 
now chosen a distance \ / 2 apart. From the same arguments as in the 
previous paragraph, the two signals are in phase when scattered back 
to the radar and add, yielding a net return signal. It is now not possible 
to choose other identical planes in the medium that have the same densily 
to cancel the signal already returned to the radar from the planes Pi and 
P2. Therefore signals are scattered from a medium in which density 
fluctuations exist in contrast to a fluctuationless medium. As a specific 
example, air contains density fluctuations. The sl^ is blue because 
light is scattered from the density fluctuations in the atmosphere. It 
is not blue when transmitted through water, because densily flue tuations 
in this medium are many orders of magnitude smaller than the density 
fluctuations in air. 

Care should be taken when considering arguments like the preceeding 
one. Although the argument is a useful heuristic tool, it is not a rigorous 
mathematical proof. One difficidiy presents itself immediately since 
planes of atoms do not scatter signals since they are of negligible thickness. 
In the above analysis, one must use a volvime in the medium which has a 
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large extent perpendicular to the incident wave and is thin relative to the 
wavelength. Then all the atoms in this particular volume of Interest 
scatter waves back to the radar receiver in a relatively coherent fashion 


analogous to the thin planes considered earlier. The two planes a distance 
X /2 apart would then be considered the boundaries of a scattering volume. 
This is, in fac t, the smallest volume which can backscatter a radar wave 
of length, \. The density fluctuations within a scattering volume or the 
density fluctuations which have scale sizes X /2 are those which scatter 
radar waves. This statement is not rigorously true if the thin planes 
being considered are thin relative to the Debye length. However, in the 
simulated case, the radar wavelength is long ty comparison to the Debye 
length and a thin plane relative to the radar wavelength does not have to 
be small by comparison to the Debye length so that it is possible to select 
planes in the medium which scatter the radar waves coherently. , 

It can be shown that the amplitude of a wave scattered from one 
particular scattering volume in the ionosphere is a Gaussian variable 
with zero mean. To prove tliis, it is necessary to consider the density 
fluctuations within the scattering volume, which generally contains a 
large number of ions. In the realm of the computer simulation, a typical 
number of particles within one scattering volume might be on the order 

g 

of 10 . The density fluctuations in a scattering volume can be considered 

g 

to be due to the random mo tions of 10 particles. The central limit 
theorem of statistics implies that these fluctuations are Gaussian about 
the mean density unless they are of such enormous magnitude that they 
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approach tlie actual average density of the medium. In this case, it 
would be impossible for tlie distribution to have Gaussian tails since a 
negative density is impossible. Given these considerations, 

A - k vn (1) 

where A is the amplitude of the signal returning from a scattering volume 
under consideration and Vn is the density fluctuation within that volume. 
Since A is a scattering wave amplitude and not a power, density waves 
within the scattering volume scatter incident radar waves coherently 
as can be seen from the arguments considering the plaiies Pi and P2. 

Since the volume scattering the wave has a particular size and position, 
the scattered wave considered in equation 1 has a specific wavelength 
and phase. This phase depends on the distance of the volume from the 
radar and upon the phase reference used a cosine wave or a sine wave. 

To understand this scattering further one can consider the phase of waves 
returning to the radar from scattering volumes. Since in the time sta- 
tionary ionosphere no particular distance from the radar is preferred 
for scattering, the phase angle of waves scattered at a particular frequency 
is a uniformly dis tr ibute d random variable . If the re turn signal s ( t ) is 
Written'' 

s(t) = A cos u)t + B sin u) t (2 ) 

for any particular frequency lu, the amplitudes A and B are Gaussian with 
zero mean since the density fluctuation in the scattering volume is 
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Gaussian with zero mean. Also, from the argument of uniform phase, 

A and B must have the same mean square. In other words, the power in 
the cosine component and the power in the sine component are equal on 
the average if no preferred phase exists in the return signal spectrum, 

It Can also be stated that the A 's and B's for each frequency component 
of the scattered wave are independent of each other in a statistical sense. 
These statements about the ionosphere may be summarized as follows: 

<AB) 0 (3) 

Consider again the limitations to equations 1, 2 and 3. One such 
limitation exists if the density fluctuation is on the order of the actual 
density in the scattering volume . It was stated that ttis is not an impor- 
tant restriction since a typical experimental radar does not receive 
signals from just one elementary scattering volume. Instead the radar 
averages over a relatively large portion of space which might contain 

0 g 

on the order of 10 to 10 elementary scattering volume s . Even if each 
one of these volumes was itself coherent, i. e. , a hard target, this 
would s imply affe c t the magnitude of the signal scattered from the volume 
of the ionosphere which the radar was examining. The hard targets 
would Ml the volume that the radar was observing in much the same 
manner that the small particles fill the much smaller scattering volume . 
Thus, for most kinds of coherent scatter, as well as incoherent scatter 
radars, the simulation being used in the program described here is 
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applicable. If the coherent target is as large as the volume of space 
being observed by the radar, the Gaussianess of the amplitudes and ran- 
domness of phases assumed in the simulation would no longer be valid. 

Thus far, no assumptions on the distribution of particle velocities within 
this ionosphere have been made. In the actual computer program it is 
assumed that the ionosphere is composed of drifting MaxY/ellians in which 
the electron and ion temperatures are allowed to differ. This does not 
affect the statistical distributions used in the computer but the fact that 
the ionosphere is of a Maxwellian character does affect the time average 
values of parameters of the spectrum. These time average values can 
be computed from theoxy once the particle distribution in the ionosphere is 
known. Given the known average values for the ionosphere parameters, 
the computer can simulate the fluctuations which must occur to 
produee these average values if the ionosphere is iricoherent within a 
radar observation volume. 

Now one can consider in more detail how the computer program models 
a signal scattered from a time -s tationary uniform ionosphere . Figure la 
describes a particular input spectrum that is fed into the computer. This 
input spectrum is caleulated from theory for the lonOsphere deseribed 
above for the following input parameters: electron temperature, 
ion temperature, T.j bulk velocily, v. ar4 density, n. These para- 
meters de termine the shape, width, and center of grayify' of the spectrum 
shown in Figure la. The spectrum is computed under the assumption 
that a radar pnils^ transmitted for ihifihite time into the ionosphere - ^ 
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and received for an infinite time try a radar receiver. Since this does 
not correspond to physical reality, the spectrum must be modified in 
the computer to look like that of Figure lb, which describes a spectrum 
that Is spread because the radar pulse is finite in length. This pulse 
is not monochromatic and the spectrum which results when the pulse is 
scattered from the ionosphere is consequently broader than that which 
would result in the ideal case when only one frequency is present. In 
Figure Ic two spectra are shown. The first is an average spectrum 
of signal plus noise power. The second is an average noise spectrum. 

This average noise spectrum is computed from 

average noise power " average signal power/snr (4) 

where snr = signal/noise ratio. Equation 4 refers to the total power in 
the signal and the total power in the noise. In order to create the noise 
spectrum of Figure Ic it is only necessary to assume that the total power 
is distributed equally among the frequency windows in the noise spectrum. 
This produces a flat spectrum which is consistent with the assumption 
of an ideal radar system. The signal plus noise spectrum is created by 
raising the signal on a pedestal so that the total power entering the receiver 
in any gate where signal is expected is the power of the signal plus the 
power of the noise. 

The spectra that have been considered thus far are smooth, averaged 
spectra. An actual example of a single spectrum (Figure Id) is selected 
at random in the computer by imposing statistical fluctuations on the 
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averaged spectrum of Figure Ic, to account for the previously discussed 
fluctuations of the ionosphere. In Figure le, the program takes account 
of the fact that the signal is received for a finite time. In the case of 
the particular simulation being considered here, the time of reception 
is equal to the time of transmission. This reception for a finite time 
causes spreading of the same nature as the spreading in Figure lb that 
is caused by finite transmission time. In Figure If, the noise spectrum 
which entered the radar uyetem from a portion of the ionosphere where 
no signal was expected, is subtracted from the signal plus noise spectrum 
which entered the radar system from a part of the ionosphere which i.<5 
to be analyzed. >yliat remains is a signal spectrum that represents 
oite possible result of the physical experiment of transmitting a signal 
into the ionosphere for finite time and receiving the signal baekscattered 
from the ionosphere for an equal finite time . This possible resulting 
spectrum can then be analyzed to determine the apparent parameters of 
the ionosphere which existed when the signal was scattered from it. 

It should be noted that the diagrams of Figure 1 are schematic in 
nature and they do not explain two important points. First, these diagrams 
do not indicate in what way spreading of the signal is accomplished in 
the computer. Second, they do not explain how statistical principles are 
applied to produce tlie fluctuations in the spectrum of Figure Id. 

To clarify these points, one must digress to consider the autocorre - 
lation func tion of the signal voltage s{t), which is defined as 
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T -T 

r ir) = TpJ s(fc) s(fc + t) dt 

■: '0 

where T is the sampling interval and t is the lag time. To illustrate 
the properties of an autocorrelation function further, consider the Specific 
function, 

s (t) = A cos (o)t + cp) . (6) 

For this s (t), 

2 T-t ■ ■ 2. ■■ 

r(T) - ^ r cos (2 tut +<UT + 2cp) dt (1 - ) cos o)T (7) 

As T equation 7 becomes 

. ■ 2, , 

r (t) = — cos tuT (8) 

■ “ ■ . .. a,, ^ 

In the finite limit when it is assumed that the phase of s(t) is uniform 
and that therefore the phase term in equation 7 averages to zero one 
obtains 

. 2 

r(T) = (1 - ™) cos CUT (9) 

Equations 8 and 9 yield ■■ ■ ' ; 

r(T) = (1-^) r (t) (10) 

This equation states that if the autocorrelation fimction is Imown in the 
limit where s(t) exists forever, it may be obtained for a signal which 
exists over a finite sample interval by multiplying try a factor (1 - t /T)* 

This multiplication factor looks like a ramp when graphed and the process 


of multiplying by it is often called triangulation for this reason. 

Triangulation is useful because the autocorrelation function '•nd a 
signal spectrum are Fourier transforms of each other. To investigate 
tills relationship, consider a signal of the form 


s(t) = Tj (a cos u) t + b sin to t) 
' n n n n 

n - 1 


(il) 


where 


(0 ~ 2Trn/T 

n 


(12) 


and T is the length of the sampling interval. Consider a lag product of 
the form 


CO CO 


s(t, )s(t_) = S Z) a a cos 10 L cos to t, 
1 ct f\ n i- n X n 4^ 


m - 0 n = 0 


H' a b COS 10 t. sin 10 t_ + b a sin to t, cos to t_ 
. m n m 1 n 2 ra n m 1 n 2 


+ b b sin 10 L sin to t„ 
m n m 1 n 2J 


(a b ) 
m n 

= 0 for all 

<a . a > 

= <b b > 

m n 

m n 

(a ) 
n 



(14) 

(15) 
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Then 

CO 

<s(fc ) s{t )> = B <a cosu) T = r(r) (X7) 

12 n n ^ 

n = 0 

where 

T = fcg “ \ (18) 

Equation 17 is the Pourier series equivalent for a real symmetric 
function of the statement that the autocorrelation function is the Fourier 
transform of the power spectrum. Notice that both the spectrum and the 
autocorrelation function are averaged. That is, the average value of 
the autocorrelation function is the Fourier transform of the average value 
of the spectrum. This is realized because of equations 14, 15, and 16, 
which are the generalization of equation 3 . If equations 14, 15, and 16 
do not hold, then the autocorrelation function of equation 17 is not the 
Fourier transform of the power spectrum, and furthermore, equation l7 
has an explicit time depc-iidence. This can be seen if the algebra leading 
to equation 17 is performed vvithout assuming equations 14, 15, and 16. 

If the autocorrelation function is an explicit function of time, that is 
depends on something other than the time difference t_ - t^ , then the pro- ^ 
cess being considered cannot be a time stationaxy process. By its nature 
a time stationary process must be independent of the time at Which one 
begins to measure physical phenomena associated with it. Another way 
of stating the same result is that if the autocorrelation function is a function 
of time, then the par am e ters of the ionosphere must be changing during 
the time that experiments are being performed. Such a situation is 


-14- 


confcrary to bhe assumptions made earlier about the characteristics of 
the ionosphere. 

To summarize, if equations 14, 15, and 16 are valid, the autocorre- 
lation function is not a function of time and is the Fourier transform of 
the power spectrum. If on the other hand equations 14, 15, and 16 are 
invalid, the autocorrelation function is a function of time, is not neces- 
sarily the Fourier transform of the spectrum, and in generai the process 
being described cannot be considered time stationary . Thus, it has been 
proved that for a time stationary process equations 14, 15, and 16 must 
hold and the autocorrelation function is a useful function to describe the 
spectral Fourier transform. 

From the above analyses, it becomes possible to spread a signal in 
the following manner. If one has the spectrum of a signal which exists 
for an infinite time and which should be spread because it will be sampled 
for only a finite time, one first derives the autocorrelation function of 
the signal by Fourier transformation. Next, the Fourier transform is 
triangulated. Then the triangulated autoeorrelation function is retrans- 
formed into frequency space to produce a spectrum equivalent to that 
which would have been derived from a signal existing for a finite time . 
The details of this equivaienee are considered in the appendix and an 
example of a spectrum before and after spreading in this manner is given 
'in.;Figime ■..2... : 

Now that the methods used for s pec tr al s pr e ading in .the o om pu te r 
have been discussed, it is necessary to understand exactly how the 
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computer produces spectra with the proper statistics to simulate the 
signal scattered from the ionosphere. Such a spectrum, estimated from 
the ISI samples obtained of the waveform s(t), will contain IS! / 2 points be- 
cause the Fourier transform of these N points produces N/ 2 amplitudes 
of cosine terms and W/2 amplitudes of sineterms. The power at any fre- 
quency is the sum of the squares of the amplitudes of the cosine and sine 
terms at that frequency. Hence, the spectral power in the ith frequency 
window is proportional to "f" B. ) and is independent of the pov/er in 
ariy other frequency window since the A's and B's are Independent of each 
other for different frequencies. 

Statistics are generated in the computer by use of 
n ■ „ : . 

u - D a ^ (19) 

where the ajs are Gaussian variables with zero mean and are independent 
of each othet'. The quantity u^ of equation 19 is a chi square variable 
with n degr{»es of freedom and is used as follows. In any spectrum, such 
as that of Figure Ic, each spectral point, whether signal plus noise or 
noise, has a sine and a cosine component. These components are squared, 
added and averaged to form the power indicated. Since the A's and B's 
are independent GaUssian variables, each frequency point is the average 
of a chi-square variable with 2 degrees of freedom . Therefore, to select 
a particular possible spectriim from the ensemble of spectra generated^ ^ ^ 

from the sea.ttering of signals from the iohosphere, it is only necesSaiy ^ - 
to select each frequency point in Fi^re Ic from a distribution hamng the^ 
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correct average value to yield the power in that point and the statistics 
associated with a chi-square variable with 2 degrees of freedom. Such 
a selection is consistent with the physical case in which one pulse is 
transmitted into the ionosphere for a finite time and received for an in- 
finite time. If one wished to generate statistics consistent with a ease in 
which many pulses were sent into the ionosphere, received> and the data 
from these pulses were averaged, one could assign more degrees of 
freedom to the distribution for each spectral point. In particular, if P 
pulses are transmitted into the ionosphere and their return spectra 
averaged, then each independent frequency point would have 2 P degrees 
of freedom. This is because each pulse yields data which is independent 
of tlae other pulses and, therefore, the number of independent variables 
being summed with P pulses are averaged is equal to 2P. 

In summaiy, then, the average spectra in Figure Ic are used to 
generate spectra with the correct statistical behavior in the following way . 
First, the number of pulses transmitted into tlie ionosphere must 
known. Second, each real frequency point is eonsidered to be independent 
of e very other real frequency point. Third, for each real frequency point 
a random number is selected from a chi-square distribution with 2P 
degrees of freedom and possessing an average value equal to ti e average 
power in that frequency window. This then generates the spectrum which 
has a s tatisticai behavior consistent with that of signals being scattered 
from the uniform, tim ionosphere already described. 
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PART m 

Data Presentation 

Two estimators are evaluated below. The first is an estimator of 
the ionospheric density and tlie second is an estimator of the bulk velocity 
of the ionosphere, which is equivalent to an estimation of the Doppler 
shift of the backscattered spectrum. These quantities are analyzed since 
they are examples of one relatively easy and one relatively difficult 
parameter to measure with reasonable accuracy using a Space Shuttle 
radar. Before presenting data from the computer simulation, it is neces- 
sary to understand how an actual error bar is computed from generation 
of many spectra such as that described in the previous section. Each 
such spectrum Is analyzed by use of the estimator being studied with 
each estimate made by the estimator pf interest being stored in the com- 
puter. The result of 1000 such estimates for a particular set of input 
parameters is presented in the histogram of Figurfe 4. In^ the simulation 
which produced this curve, a spectrum with the parameters and the 
signal-to-noise ratio listed in the figure was given to the computer as 
input data. The computer then produced the given histogram having an 
apparently Gaussian shape and a maximum at a mean value of almost 
200 meters /second. The actual mean value estimated by the computer 
is 196 meters /second and the standard deviation is 32 meters/second. 

Now that the process producing the standard deviation or error bar 
has been described, the simulation estimates of errors in density and line 
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of sight velocity will be discussed. It is possible to develop a theoretical 
equation for the fractional error of the density, An/n in the following way. 

where P, is the total power, the fractional error in the density may be 

t ■ " " ' ' 

obtained from the fractional error in the return signal power. To find 
the fractional error in the total power consider one point in frequency 
space, whose power is p. If the point contains signal and noise power then 


( 20 ) 


(p ) = c (1 + — ^) 

*^s + n snr 


( 21 ) 


where snr is the signal to noise ratio and c is a constant that represents 
the average signal power. Hence 


<P.) = 


n snr 




( 22 ) 


Equations 21 and 22 describe every point of a real signal plus noise and 
noise spectrum if both spectra are considered to be flat. Both the signal 
plus noise point and the noise point are chi-square variables with 2 degrees 
of freedom. For this type of distribution it can be shown that 


, 2 , 2 . 

Au = <u ) = 


( 23 ) 


.2 


where u is a chi-square variable with 2 degrees of freedom, and Au 
is the variance. Since the variances of independent quantities are additive 


2 2 2 
A (u^ - Ug) - A u- 


if u, and u„ are independent variables. 

' 1 . £t 

In order to obtain a point in the final signal specti’um as in Figure If 
it is necessary to take a point in the signal plus noise spectrum and sub- 
tract from it a point in the noise spectrum. From equations 21 through 24 


A (P . - P [ (1 -I- — )^ 4 (— )^ ] 

s + n n *• snr snr r 


(25) 


From equations (25) and (22) 
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where is the final power in the particular point being considered. 

Equation 26 applies to the power in only one point. However, there 
are N/2 independent frequency points. Each one of these points contri- 
butes to the power whose estimate is therefore improved by a factor of 
I//N/2'. In addition, whenP pulses are transmitted and the data from 
the different pulses are averaged to yield a spectrum, another improve - 
ment ty a factor l//^ is accomplished. Thus, the final form for the 
fractional error in signal power, or from equation (20), for the fractiohal 
error in the density estimate is 


# " t(i + (I® / 21 


n 


snr 


snr 


Since this formula was derived under the assumption that the signal spec- 
trum is flat and since signal spec tra of interest are hot fiat> the actual 
errors are expected to be somewhat greater than those predicted by this 


theoretical formula. This is because all of the points do not contribute 
equally to the variance of the power. To give an extreme example, if 
all of the power in a signal were concentrated in one point, then there 
would be no enhancement of the accuracy by a factor 1//N/2' as is assumed 
in e qua tion 27. Howe ve r, for mo s t s pec tra, the d e via tion of the signal 
spectrum from flatness is expected to increase errors by a factor of 
only about 10%. Thus, the theoretical formula for fi’actional density 
error derived in equation 27 is expected to be accurate over wide ranges 
of signal to noise ratio and numbers of pulses. 

To verify this statement, consider Figure 5 which is a plot of An /n 
estimated for a range of ionospheric parameters by the computer simu- 
lation versus that given by equation 27. The dots of this figure represent 
a high temperature case having an electroxi temperature of 3500° and 
an ion temperature of 2500° while the crosses represent a low temper- 
ature case with equal electron and ion temperatures of 1000° , If a case 
where the signal to noise ratio and number of pulses is high Were sinau- 
lated'in the computer, this produeed a point lying near the origin corres- 
ponding to small values of An/n. The pdints representing small signal 
to noise ratios and small numbers of pulses lie in the upper right hand 
portion of the graph. 

If the equation agreed perfectly with the simulation, every dot and 
eveiy cross would lie on a line a t an inclination of 15 ° . As can be seen 
from the figure, the dots do lie close to this line over tlie wide range of 
parameters used in the simulation. However, it should be noted that the 
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crosses lie consistently above the line, indicating that the theoretic. * 
equation consistently underestimates the density error associated with 
them by a typical factor of about 10%. This result is explained by the 
narrowness of the low temperature spectra associated with the crosses. 

In summary, the computer simulation for broad spectra show that 
equation 27 is an excellent predictor of the fractional error in the density 
estimate. For relatively narrow spectra with low electron and ion tem- 
peratures, it underestimates density errors hy about 10%. Since the 
feasibility of a space shuttle radar does not depend on 10% corrections 
to such error estimates, equation 27 will be used to evaluate the perfor- 
mance of such a radar. 

In Figure 6, contours of constant An/ n as determined by equation 27 
are plotted as functions of the signal to noise ratio (ranging from 0. 1 
to 100) and number of radar pulses, P (ranging from 1 to 10 ). 

Again this graph is a log-log plot as in Figure 5. A general char- 
acteristic of the contours in this figure is their tendency to parallel 
the snr axis for signal- to -noise ratios greater than about 2. This 
break point at snr « 2 implies that the accuracy of a densily measure 
ment is not greatly improved by making the signal -to-noise ratio much 
greater than 2, 

The break points in the contours can be understood by considering 
that beyond some value, no further increase in the signal- to -noise ratio 
will produce a me aningful decrease in the number of pulses for the same 
value of An/n because there is error associated with the signal power 




eve;H in the absence of noise. When a value of the signal- to- noise is 
reached where most of the received power is signal, most of the error 
associated with the measurement is due to the statistical fluctuations in 
the return signal and is unrelated to the noise. These statistical fluctu- 
ations exist because the signal is scattered from a alatistically fluctuating 
medium. 

It is thus concluded that there is little point in designing a radar- 
system having a large signal- to-noise ratio. After an optimum signal- 
to-noise ratio is reached, one can decrease the error bars in a measure- 
ment only by increasing the number of independent samples taken. For 
purposes of discussion, an optimum signal/noise ratio is considered to 
be- 2 . 

The analysis presented above yields some useful order of magnitude 
numbers. The density can be measured to an accuracy of about 10% 
with about 10 radar pulses if the signal- to -noise ratio is bigger than 
about 2. Density measurements to an accuracy of 1% require about 1 000 
pulses. If the radar operates at 100 pulses/second and the space shuttle 
moves through the ionosphere at a speed of 1 0 km/ sec, a measurement 
requiring 10 pulses would occur over a distance of 1 km. A measurement 
requiring 1000 pulses would involve 100 km of spatial averaging, which 
might destroy the utility of some types of measurements. 

The velocity estimatox’ will be examined by developing an equation 
for the uncertainty in the velocity analogous to equation 27. Since, at 
least in some range, the velocity estimator might be linear in the density 
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estitnator, an attempt was made' to design an empirical equation which 
would fit the velocity errors as closely as possible, remains simple, 
and have a form relatively similar to that used for the density estimator. 
The equation which resulted from this process is 


i.Wsec) = [l+(|)“][(l+5^) .(-.y ] 


(28) 


This equation was developed through analyses of data such as tliat 
presented in tables 9 and 10. Table 9 compares the signal- to -noise 
dependence of equation 28 with the actual dependence of the average 
data. The columns in this table are veloci^ errors estimated from the 
computer simulation, for different pulse numbers and two different temper- 
ature cases. The rows in the table represent different signal- to-noise 
ratios. The last column is an average of all the numbers in a particular 
row, with values in excess of 475 m/sec dropped from the table. This 
was done because the computer program has a cutoff Which will not allow 
it to estimate any velocity greater than 800 m/sec. This means tliat 
velocity errors which approach this value are , progressively underestimated 
and these data were eliminated from the table. The table values are 
normalized to the error associated with a signal-to-noise ratio of 100. 
Comparisoh of the average signal-to-noise dependence with that predicted 
by equation 28 produees close agreement, even though there are uncor- 
rected temperature dependences in the data of Table 9. ^ 




In Table 10 the pulse dependence In qquation 28 is compared with that 

of the computer simulation. Here the error is plotted in units of 

Av,,/Av„„ where is the velocity uncertainty for P pulses. Ex- 
P 22, 500 P 

amination of this table indicates that the equation predicts the pulse 
dependence for large numbers of pulses and overestimates it for small 
number of pulses. By examination of the tables the statement that 
equation 28 predicts velocily errors to within 15% typically and in the 
worst case to 30% accuracy can be verified. 

Figure 7 presents a graphical comparison of equation 28 and the 
velocity uncertainties estimated by the computer simulation. Koto that 
the points do not lie as close to the 45 degree fit line as they did for the 
densily errors of Figure 5. Further note that not all of the dots and 
crosses lie consistently on the line and that for example sometimes 
cx’osses lie below the line while other times they He above the line. This 
shows that in the estimation of velocity there are temperature dependences 
which are not bemg considered here. Note, however, that the overall 
accuracy of the equation is (ypically 15 or 20%. Since 15 or 20% errors 
in the uneertainly of the veloetty estimate are small enough to be neglected 
when obtaining rules of thumb for operation of the space shuttle radai-, 
such ei’rors will be neglected. 

Contours of constant velocity uncertainty are presented in Figure 8 
as functions of the signal- to-noise ratio and number of radar pulses, 

These contours show many of the features of the density equation cobt-gfars 
since the equations which produced both sets of contours are relatively 
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similar. In both curves, the break point for improvement of accuracy 
versus signal- to -noise occurs at a value of about 2 for all contours. 

For signal- to -noise ratios in excess of this value, one can obtain an 
accuracy of 100 m/sec in the velocity estimate with appi^oximately 15 
pulses, wMch involve a spatial averaging of about 1. 6 km. If one wished 
to obtain an accuracy of 10 m/sec, one would need to collect approximately 
1500 pulses. Tills, would indicate spatial averaging over an area of per- 
haps 150 km, which may be too large for many classes of experiment. 
Thus, for experiments in the auroral zone, where typical bulk velocities 
of tlie ionosphere are much larger than 100 m/sec, tlie space shuttle 
radar would be a useable tool for analysis. In the equatorial zone, how- 
ever, where ty’pical velocities are less than 100 m/sec, the error in the 
measurement would be comparable to the magnitude of the velocity being 
measured unless one averaged over a large volume of space. Therefore, 
Space shuttle experiments in the equatorial zone designed to measure 
bulk velocities of the ionosphere are probably not feasible. 

In summarzing the data the following conclusions can be made. Given 
reasonable pulse repetition rates (perhaps 100/sec) a sufficient signal- 
to-noise ratio (a 2), aiid a typical space shuttle speed (10 km /sec), one 
can make density measurements to the accuracy of 10% and velocity 
measurements to the accuracy of 100 m/sec, over small spatial intervals. 

Certain cautions should be kept in mind when extrapolating from the 
data presented in part III. It must be remembered that the computer 
simulated a particular operating mode of the Chatanika radar. In this 
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mode the radar transmitted a pulse for 320 jisec and received a pulse 
for the same length of time. This corresponds to a range resoJution of 
approximately 50 km. During the received pulse, 32 samples were taken 
= 32, This results in a bandwidth of 50 kHz containing 16 independent 
discrete frequency points. When one extrapolates from the data in part 
III these explicit dependences should be kept in mind if one wishes to do 
analyses for different range resolutions or different sampling rates. 


APPENDIX 


Throughout this appendix the following formalism will be followed. 
s(t) will represent a continuous function, s^. will represent a discrete 
function existing only for specific values of t. It will also be assumed 
that the signal under consideration contains no d. c. nor Nyquist term, 

^0 ~ ^N/’ ~ terms are neglected because they complicate the 

mathematical formalism and add nothing to the understanding of discrete 
vs. continuous analysis. 


I. Discrete Functions and Qrthogonalily 

Beginning with continuous functions, one can study how the equations 
of Fourier analysis are modifisd as discrete sampling is Introduced, 
if on an interval T 

s(t) = S a cos ncD^t + b sin A (1.1) 

■ n ■= 0 ' ; , : ; 

w here u) ^ “ 2 rr / T , then for all n 


y-' T ■■■ : 

2 

a = 7 = r s{t) cos mu-t dt 
n T ‘’p 0 


( 1 . 2 ) 


and 


- I" f s(t) sin nu),,t dt ^ (1. 3) 

n ■■ -'T ■ : ■; 0 

These equations follow directly from equation 1 because of the ortho- 
gonaliiy of sin and cos on the interval T- To form a discreteisignal related 
to s (t) one can divide the interval T into N in ter vals of length T /N , with 
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the first interval labeled as the 0 interval and the last interval labeled 
as N - 1. Imagine that at the beginning of each interval s(t) is sampled 
and its instantaneous value stored, This converts s(t) into a new function 
S|, with N discrete values on the sample interval. This new function has 
the value 0 except where 

k = 0, 1, , . . . , K -1- Sj. can be described mathematically as 

■ : T ■ 

^ (1.5) 

'■or " 

k “ 0 , 1 , 2 , « • • j N “■ 1 . 

Once the interval length T and number of samples N are defined, 
equation 1. 6 is a function of k alone. For this reason, the discrete 
signal is often called s^^. 

At this point it is natural to inquire whether a signal of the fype in ^ 
equatiGii 1. 6 can be represented in a manner analbgous to the represen- 
tation of s(t) in equation 1.1, In fact, such a representation is possible 
but the pr oof of this point will be pursued indirectly . Before the proof 
can be imderstoodr one must inq^ of equations 1, 2 

and 1.3. TheseequationSdependfortheirvaliditybothQntherepre-*i \, 

sentability of s(t) as in equation 1, and on the orthogocailiy relations 
for sin and cos terms in Fourier series. 
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The orthogonality relations in the continuous case are 


if 

V J 


T 


T 


cos mwQt cos nu)Q 


t dt 


1 ^ 


Sin moj. 


t sin nouQtdt 



m, n 



T 

Tjr r COS mw- I sin n<o„ tdt = 0 
T 0 0 


(1.7) 


for ail m and n. 


The orthogonality relations of equation 1 . 7 can be derived from 


the trigonometric identities and die fact that ary sin or cos wave periodic 
on the interval T has integral 0 unless the term being considered is a 
constant, that is, unless the term being considered is a cos term with 
argument 0. As an illustration consider 

rj> 

^ J sin mOJQ t sin nuj^ t dt = 


T - :: , ■ - 

J cos (m - n) (U -t dt - 

n . ' V 


T 

J cos 
0 



1 A 

2 A; 


m, n 


( 1 , 8 ) 


where m and n are positive integers. 

If in the discrete case one were to discover that periodic functions 
summed over an interval behaved similarly to the periodic continuous 
functions one could prove the orthogonality relations at once, 

If One considers an interval of length T divided into N parts then in 
acGordanee with equatibn 1. 6 the orthogonaltty reiationS in equation 1 . 7 
for discrete s^, become 



™ S [cos mu)Q (k T/N) cos nWj, 
■ k = 0 '■ 


(k T/N)]At = 



B [siji mu>_(k T/I^) sin nu)_ (k T/lNi)lAt ^ 
k-O*- ° ° 


t A ■ ■■■ , ■ 

TO -.n, pN 

' r'^'K . ■' ' 

B [oos ma)Q {k T/'N) sin nojQ (k T/N)jAt = 0 (1.9) 

" 'k;" O' ' 

Here, A t equal T/I^ is the length o^ S- sarople interval and is analogous 
to the diTferential dt in the integrals of 1. 7* From the definitions of 
At and u>g equation l> 9 can be written as 

;,""N -'I' 

rr B cos{2nmk/N)cos(2TTnk/lS}) = 

V" H-i": 

— B sin(2nmk/N)sioX2unk/N) = 

■'■k'=^o'" 

^ B cOS(2tt nik/]SI)sin{2TTnk/N) = 0 (1.10) 

■r;''" V , k'=0':x^" X, / ■; 

where p is ah integer . 


To derive the orthogonality relations for discrete signals one needs 
only trigonoTOe trie identities and complete knowledge off and g defined as 
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N-1 

X = -rr Tl COS 2 tt qk/I^ , 
9 N ,,^n 


N-1 

ff = — S sin 2n qk/K 
^ ■ k = 0 


To find f and g consider the relation 

q q 


1 \ 
f^ + i = :jT ( S COS 2Trqk/N + i sin 2rr qk/N j 

^ ^ ^ A k = 0 / 


1 i2TTqk/N 

ISJ ^ ® 

k = 0 


Rewriting of equation 1. 12 yields 


N -1 


V‘«,-w A 


k = 0 


By use of the algebraic identity 


k I-x^’ 

X - —— — 

k = 0 


Equation 1 .13 becomes 


r * _ 1 1 - V e J 

q '■^q H ^ _ ^i2TTq/lS! 


q ’■^q N . atriq/lSI 

V. . . . , ■■ . : 1 - e 


This gives 
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f +ig = 0 

q ®q 


q f pN p = 1,2, . . 


If q = pN, then both the numerator and denominator in equation 1.16 are 
0, and one must use equation 1.12 with q = pN. 


W-i 


^ is ~ 77 ^ cos 2rrpk + i sin 2TTpk = 1 


’q N 


k = 0 


where q = pK and p and k are integerSi Equating real and imaginary parts 
and combining the results of equations 1.17 and 1.18 yields 

trr B cos 2TTqk/K =1 

■k=-0:.\ , , ■ 


q = pN - 0 otherwise. 


N- 1.. : . 

B sin 2Trqk/W = 0 for all integer q- 

-k=:0' ' VV; 


These equations are analogous to the integral relations 

,.T ■ ■ T 

— J cos moj- t dt = J cos dt = 1 


if m = 0 and 0 otherwise 


~ J sin miu t dt = J* sin — dt = 0 for all m 


Wote here the periodicity in equation 1 il 9 which does not exist in equation 
1. 20. This leads to periodic behavior in the orthogonality relations 


equations 1. 10 for discrete functions which does not e#s t in the case of 
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Fourier series intended to analyze continuous functions. This problem 
of periodicity will be discussed in the next section along with the problem 
of representabiliiy of a discrete function by a Fourier series. 

II. Periodicity and Completeness in Frequehcy Space 

This section contains explicit expressions for a and b for discrete 

■ ■ n n ■ 

functions analogous to equation 1.2 and 1.3. These expressions are derived 
using equation 1 . 10. However, because of the peculiar periodicity exhibited 
by equation 1.10, confusion about these results may occur unless efforts 
are made to understand in what way the relations for a and b are 
periodic. 

if one writes expressions for discrete functions a^ and b^ naively, 
one obtains 

„ N-l N-1 

: .. n-t.., , ; 

n IS! Z) s, sin 2Trnk/I\' (2.2) 

Here"..':'.': 

-gU'/T : 

• ■ ^ At/='T/]SI: / " ; ; \ ^ ;f2. 3) . 

where t is the length of the sampliug interval, and N is the number of 
samples taken in the interval. Relations 2.1 and 2. 2 give 
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3 . 3 . , 

n n+pW 


b = b + pN 
n n ‘^ 

p = 0,1,2... 

and 

a = a „ 

n pIS! - n 

b = -b 

n pN “ n 


(2.4) 


(2.5) 


Setting p = 1 gives 


a 

n 


W - n 


b 

n 




( 2 . 6 ) 


It follows from equation 2. 4 that if Uq. . .a^^ _ ^ and bg. . . .bj^ _ are 

Icnown, then a and b are known for all n. It follows from equation 2. 6 

, n n .. 

that if a and b are known then a^, and b^, are known. Thus, one 

:n ■.;,n'.. , •■■N-'- n' N.-.n: ;. 

need know only a^. . . a^ and bg. . .b^ where r is the greatest integer less 
than or equal to N / 2 in order to completely determine the frequency com- 
ponents a and b of s, for all n. 

Writing an expression analogous to equation 1. 1 now yields 
00 

s, = S a cos nou- (kT/W) +b sin niu- (kT/N) = 
k ' n 0 n 0 

n = 0 



2 TT nk . , . 

2j a cos —irr— + b s in 

_ Q 


2frnfc 

: w ■ 


(2.7) 
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One may well ask at this point if it is necessary to extend the sum 

in equation 2. 7 from 0 to co in order to completely represent Sj^. It is 

already known that the a's and b’s repeat beyond the index r. If the functions 

in equation 2 . 7 multiplying the a^'s and b^'s repeat as do the coefficients 

a^ and b^, then additional terms in the sum beyond r in equation 2. 7 

will contribute no new informabout about the functional form of s, . To 

decide whether the ?um in 2.. 7 is periodic, consider one term in that sum, 

the mth term, where m is between 0 and r. Gall the contribution to the 

function s, of this term, s, . Then 
Ic K* m 


2Tr(m + phi) k , 

s, . TV 1 “ S' j TM cos rr + 

k, m + pN m + pN N 


, . 2TT(m + pN) k 2rrmk , 

b / sin — rr— - — = a cos — — + 

m + pN K m IS! 


, . 2rrmk 

b sin — r; — = s, 

m N k, m 


Thus result follows from equation 2. 4 and uses the same techniques in 
its derivation as were used in equation 2. 4. Further, 

2n(lS! - m)k , V • 2Tr {N - m)k 

s, ,,, = a., cos r; — — +b„ sin - — 

k, IS! - m N - m N N, - m IS! 


cos (2tt -2-rrmk/lS! ) + b_, sin (2rr - 2wmk/W ) 
jS! - m In - m 


a- cos - b^. sin 2 trink/N 

IS! - m K N - m 


2TTnnk , v . 2rrmk 
a. cos '■ ■■ ■■ + b sin . 
m N m N 


k, m 


( 2 : 2 ) 
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Combining the results of equations 2.8 and 2« 9 yields 


®k,in ®k, pN * m 


Equation 2. 10 states that if s. is known for 0 < m < r then adding terms 

. ■ . ■ , ZXX 

of the form on the right hand side of equation 2. 10 to corresponding terms 
of the form on the left hand side of this equation for all values Of m 
between 0 and r simply multiplies the sum in equation 2. 7 by a constant 
but does not change the form of Sj^. Thus, without loss of generality 
one can write 

vv 2TTnk , , . 2rrnk , 

s, = S s, - Vi a cos ■ - ■ • + b sin — — ( 

k - k, n « n N n IS! 

n = 0 ' n = 0 

Generally, it follows from equation 2. 10 that the limits on equation 2.11 
can be written as 

^ 2TTnk , , 2iTnk 

s, = Li a cos —t— ~ + b sin — — = 
k T .1 u N n N 

: . , y,. m,=,plS!, , \ 

^ 2TTnk j ^ . 2nnk / 

Ij a cos— — + b sin — ^ — ( 

, . T,, n N n N 

n = r + 1 + pN 

where p is any integer. Equation 2. 12 shows that, given a particular Sj^, 
it is impossible to determine its frequency components. Specifically, 
for pv=;0 


2rrnk , , . 2Trnk 

s, = 2L/ a cos ~— + b sin — rr— 

k „ n K n N 

. n:=0 : 


w,_ r : 

^ ^ 2TTnk . , , 2Trnk 

2j a cos b sin — irr— 

. , n ^ N n ;v E ■ 

n = r + 1 
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SLnce, r s N/2, equation 2.13 states that one cannot decide whether 
has frequency components between 0 and K Wq/ 2 or between Nu>p/2 and 
WtUQ. Thus if one samples at a linear frequency l/At = l/(T/IS!) one 
can resolve components unambiguously only up to a frequency of 
Ku)q/2 X 2tt = N/2T - 1/2 At, Thus, the maximum frequency which can 
be unambiguously determined equals one half of the sampling frequency. 
This is a proof of Nyquist's theorem for discrete functions. 

If on a fixed length sample interval N -* m the number of a^'s and 
b^'s one must obtain before the coefficients repeat also approaches 
infinity. This is reasonable since in the limit N -• m Sj. "* s(t), a contin- 
uous function. In other words as the struetui‘e of Sj^ becomes finer, rapidly 
oscillating sine and Cosine terms must be added to the summation in 
equation 2.11 to represent the rapid fluctuations in the fine structure. 

A plausability argument for the fact that a finite number of terms is 
needed to represent a discrete function is that one would expect a function 
with a finite number of points to be representable with a finite number 
of functions which were independent on the interval of summation. 

Now that it has been demonstrated that Sj^ need be represented with 
only a finite number of sine and cosine terms, one can show that e very 
function s, can be represented with a series of the form of equation 2. 11, 
as follows. It has been noted that equations 1. 1 9 and 1 . 20 are similar 
but that equation 1.20 possesses a periodic behavior. If one is limited 
to a series of the form of equation 2. 11, the periodicity is removed 
from equation 1.20 and equation 1.19 and 1. 20 imply identical orthogonality 
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relations for discrete and continuous functions. Thus, equation 1.10 
becomes 


1 

N 


N - 1 

n 

k = 0 


cos 


2Trnk 

N 


cos 


2nmk 

N 



N -1 
S sin 
k=0 


2nmk . 2nnk 



A 

m,n 

(2.14) 


W - 1 

1 V* 2TTmk . 2miK 

7T 2./ cos jp- sin — 

k^O 


for 0 < m < r and 0 < n < r 


In order to decide whether or not any function s^^ can be represented 
in the form equation 2. 11, it is only necessary to represent one particular 
point s_. of that function. The question of representabilily is equivalent 
to the question of completeness and can be stated in the following terms. 

If one computes the Fourier components of a function and then sums these 
Fourier components does one obtain the original function? If the answer 
is yes for all functions, then the set of Fourier components is said to be 
complete. In the case of the discrete function, if one can represent any 
one point in the set completely, one can represent any function by summing 
the representations for individual points by the principle of superposition. 

Let s^. be any point in the sample space. Define 

S, = cA (2. 17) 

■ k k, .1 
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which states that jth point on the sample interval has a value c. 


N-1 


a . ” S C A, . 

n» J N k.D N 


2nnk 


u - 2 ^ , . 2irnk 

. = -T 1j C a, . sin • 

n, 3 ^ k = 0 ^ ^ 


_ 2c 2Trni 

a . = -rr cos 
n, 3 N N 


1 _ 2c . 2Trni 

b . = rrr sin 

n, 3 N K 


Using equation 2.11 


^ _ - 5 -v t 2c ^ 2 Trn 3 2r,nk' , ijirnj . 

s,. " 1 cos ^ c os + s in sin 


2iTnj . 2TTnk 
_ s,n _ 


Extendirig the interval of integration from 0 to N - 1 simply doubles Sj^ 
by equation 2.13. Therefore, 


K - If 

e 2TTni 2TTnk ^ . 2TTn1 . 2:Trnk 

\ ° M N + 

■ n= O'L . 


which gives 


■ N-1 * .. , . 

c 2 Tin (3 - k) 

- N : 

n = 0 


By equation 1 . 19 


s, = cA, . 
k k.3 
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Equafcion 2,23 is the equivalent of the completeness relation for Fourier 
integrals. Its existence proves that can be completely represented 
in the form equation 2. 11, 

Thus, a representation s, in terms of a discrete Fourier series has 
been found. Further, this representation was found to be periodic in 
its coefficiente. Lastly, it was proved that s^^ can be Completely represented 
by allowing the index in equation 2. 11 to range from 0 to r only. If one 
inspects the iiumber of terms thus included in the sum in equation 2. 11, 
one discovers that the number of coefficients in the Fourier series is 
equal to the number of points to be represented. This is reasonable since 
all the coefficients multiply functions which are independent on the interval 
of summation, so that IS! independent points are represented by linear 
combinations of N independent functions. 

lii. An Analysis of Spreading by Use of Fourier Series 

In sections land II of the Appendix, the mathematical formalism 
used to derive amplitude spectra for discrete signals was developed. 

IMext, certain applications of this formalism will be considered. The 
first of these is the uncertainfy principle which expresses itself as a 
spreading of the signal width in frequency space when the signal is sampled 
for a finite time. It will be assumed throughout this section that the 
power in the signal s pectrum is proportional to tlie square of the sum of 
the Fourier amplitude s for a given frequency just as is assumed for 
continuous representations using Fourier series. 
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Consider a signal of the form 


2TTiJ.k , . 2iTiJ-k 

= IT cos + V sm 


where the variable it need not be an integer. Consistent with previous 
definitions, the discrete representation of a wave with linear fre- 

quency p/T* An alternate representation is 


I 2TT|ik , .. 
Sj^ - a cos l 


^ ^ 2 2 ^ 

where cp = tan (-v/u) and a = (u + v )^. Since p is not an integer, Sj^ 
in equation 3. 1 is not periodic on the sampling interval, From equation 
2.1 and 2. 2 


■ Is! - 1 ' / 

2 2nnk / 2iTiik . . 2TTak 

a = r- 1 j cos — rr— lu cos— 7— — + V Sin — 

n N V n N \ N K 

k = 0 


N - 1 / 

2 -O • 2 link ( 2TT:|tk , . 2TT|tk 

... k = 0 ' 


Define a new quantity c - a -f ib ^ The spectral power per unit frequency 
^ n n n 

. / 2 . . 2 \ „ .. I I 2 ; « , rni 


is p = ( a b \/2u)„ - j c 1 /2uj_. Then 
^0) A n n / 0 n' 0 


2 ^ i(2Trnk/N) 

C : ' ■ 2-j e ■ 

” ® ic=o ' 


2Trpk , . 2Trpk 

u cos ^ + V sin 


I 
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c 

n 


i. 

N 


N-1 

S 

k = 0 


L(2TTnk/W) 

G 


11 ( e 


i2TTprk/W 


•h e 





^“i2TTpik/i'!^j 


. _ 2TT(n - |ji) 

^ - N 


V - 


_ 2n(p + |J.) 


ISi 


Eqtiafcion 3. 6 becomes 


N -1 


1 „ / Lvk , lAkV . / ivk iAk\ 

c = ^ B u (e V + e ) - IV (e ^ - e ) 

n ISI , V f \ ' 


k = 0 


With the algebraic identily 1. 14, equation 3. 9 can be written 




1 - e 


ivN 


n N 


1 - e 


lY 


+ (u + i^ 


\ iAW 
1 - e 

ii J 


1 - e 


earlier, assume that u and v are Gaussian variables with 

'■■■■ '-Xu)'' V) '■= 0 ^ 

/ 2. y 2» 

(u- ), - <V > , ■ 

<uv) = 

This is a mathematical statement of the randomness of the phase and 
amplitude of a signal s^ scattered from the ionosphere. Define 


(3. G) 

(3.7) 
(3. 8) 

(3.9) 

(3.10) 


{3.11) 


-A17- 



P = - — (3»12) 

■■ ■■■ 1 - e^"' 

Then equation 3.10 can be written as 


Cn = ^ [o: (u + iv) + P(xi - iv)] (3. 13) 


1 c I ^ [ 1 aj (u^ + + aP * (u^ - v^) + 

a^p (u^ - + \ ^\^J (3. 14) 


where P* is the complex conjugate of P. With the use of the physical 


.11 


/ 1 1 2v 1 / 2 , 2s / 1 - ( 

Clc^l > <U 


2v 7 1 - COS N A ^ 1 “ cos 

'T' 


cos A 1 - COSY 


(3.15) 


For the sake of explieitness this eg-uation can be written as 


1 2i 1 t 2 2» 

Xi . w 


N 


1 - cos 2Tr(n + ii) 


1 - cos 2rr(n - p,) 
1 - cos 


1 - cos 2ir(n + 


(3.1 6) 
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Equation 3.16 is proportional to the average power per unit frequency in 
the nth frequency window. Note that one term in the equation is symmetric 
about the frequency |i/T while the other term is not. This means that 
if a signal which is not commensurate with a sample period is sampled 
for a finite time, it will be spread and this spreading is consistent with 
the uncertainty principle. However., the signal will not be spread symme- 
trically about its original frequency. 

Before 3.16 can be applied to a spectrum, one must know the conditions 
under which equation 3. 16 applies to a spectrum rather than simply to a 
monochromatic signal. It is possible to imagine cases in which cross 
terms would occur between different frequencies in the average power 
spectrum of equation 3. 16 making it invalid for continuous spectra. To 
pursue this point further, assume a signal of the form 

s, = S u cos 2tTii/lS! + V sin 2ttu/N (3. 

k u, u, 

■■li ■■ - ,■-■■■■ ■ .v,./: 

The sum over p is purely formal since p, is not necessarily an integer. 

If p, is allowed to become a continuous variable, then Sj^ would correspond 
to a continuous signal and the spectrum would likewise be continuous. 
However, the subscript notation is retained for clarity and to preserve 
continuity with prior results. Under the. conditions consistent with a 
time stationary Gaussian process of the type discussed earlier. 
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(u V ) = 0 for all li and v 

Then equabion 3. 17 applies bo a specbrum and we ean wrlbe 


(3.18) 



1 ^ / 2 2, 7 1 - cos 2fT(n - |J.) 

p, ^}A ' 1^1 - cos 2-rT(n - |j,)/lSI 


1 - cos 2tt (n + \i) 

1 - cos 2n (n + p.)/1SI 


(3.19) 


This equabion remains valid because bhe condibions of equation 3. 18 insure 
that cross-terms from different jj, values will not contribute toX j c^| ) • 
Since it is now indicated that equation 3.16 applies to continuous spectra 
in the form 3. 19 when the process being considered is Gaussian and time 
stationary, one can draw certain conclusions about the spectra associated 
with such processes. The most important of these conclusions involves 
the eenter of graviiy estimator. This estimator attempts to determine 
by what amount tlie center of the spectrum is shifted from a position 
symmetrically placed abbut the center of the bandpass. There are indi- 
cations from the forms of equation 3 . 1 6 and 3 . 19 that the s pectrum after 
being sampled for a finite time might have a center of gravity shifted 
relative to the center of gravity when the spectrum existed in a continuous 
state; This Mnd of shift could lead to systematic errors in the center 
of gravity estimate. The expression 3. 19 for spectral spreading is an 
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extremely interesting one and will be discussed at much greater length 
in the next section. However, before this discussion it will be shown 
tliat the expression 3. 19 obtained by squaring of Fourier components can 
also be obtained by Fourier transformation of the triangulated autoeorre - 
la tion function, This is strong indication that expression 3.19 is a gen- 
eral characteristic of finite discrete sampling rather than an apparent 
characteristic resulting from some particular method of obtaining tlie 
spectrum. 


To show that triangulation, the process discussed for continuous 
functions in equation 9, achieves the same effect as does the squaring 
of the Fourier amplitudes of a finitely sampled discrete signal, consider 
a discretely sampled function of the form 
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N -k-1 

r = rr S s .s . , 

k N ] 3 - l- k 


Substituting equation 4. 1 into equation 4. 2 gives 


1 ( r 

S 11 cos 2TT(ij/N + V sin 2ntij/N 

i = 0 ^ ib I- ^ ^ 


Jj u cos Zrrvtj + k)/N + v sin 2 ttv (3 + k)/iS! 


‘ i ssfu u COB %ii cos 
k .^ Nm-vuV ]SI N 

D-0 Lr 


■ . 2Trv{i+k) , ,2TTaj 2Tfv(3 + k) 

+ U V cos sm^^^ ^ -^ + V U SU1-— T^ COS- — 

a V N w a ^ N N 


, . 2ttvj . 2na(] + k) 

+ V V sin ■ " - sin ■ 

a V N N 


If one imposes constraints on the u '5 and v 's consistent with those of 
equations 14, 15 and 16 


2 2 

IS! - k - 1 u + V 

<rk> = S D < ^ ^ > cos 2rrvk/]S! 


Noting that equation 4. 5 is symmetric in k and contains no j dependence 
it can be written as 


<r^> = 


I ^ I \ U + V 

)Z/ < ^ ) cos 2Trvk/N 

•JN. s) ■ ■ 
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From equation 17 it is known that r(T) is tlae li’ourier cosine transform 
of the power spectrum. It must be remembered however, that the auto- 
correlation function and the power spectral density are cosine transforms 
of each other only when they are real and symmetric. Therefore, one 
must perform tlie cosine transform over both positive and negative values 
of its argument. Otherwise, a real symmetric function would not have 
a Fourier transform which was real and symmetric. The fact that the 
transform is two-sided, that is, performed for positive and negative 
values of the argument, will become important when the discrete case 
is considered. 

The discrete relation corresponding to equation 1 7 is 


r nuj^kT 

<r(kT/lSI)) = S <p > cos — TT — £iu) 

n=0 ; • ^ ^ 


where 





Under the assumption that equal power s exis t in the sine and cosine 
components. Equation 4. 7 can be rewritten as 


r a +b 

<r > = Z) <-2 — - — S_ ) cos 2rfnk/N 
n=0 ; 


(4.7) 


(4.8) 


It should be noted tliat this function is symmetric in the variable k as it 
should be since it is the Fourier transform of a real symme trie spe c trum. 
Because (r, ) is symmetric, it contains only cosine components and can 

■ K. 

be written as 


-A23- 


<r. > = S A cos 
k 

n = 0 


Equations 4. 8 and 4. 9 are identical in form if is identified with the 
spectral power in a particular frequency window, A^ = ^ 

Thus, the spectral power can be obtained by determining the A^, Sub- 
stituting equation 4, 6 into equation 4, 9 leads to 

- 2 , 2 

A - '"v V 2nvk 2unk 

\ n S <l - lkt /R)^ cos cos — ^ 


To simplify notation rewrite equation 4. 10 as 


^ /U + V V 2 ^ , /-N.i 2Trvk _ 2rrnk ,, 

k = Q 

where the index v has been suppressed and notation is consistent with that 
used for the one component case. 

Defining delta and gamma as in equation 3. 7 and 3. 8 with jj, in place 


of V gives 


' 2 . ' 2 ' ; 

A = < - ^ o-m " ™ ') S {1 - k/bJ ) {cos Ak + cos Yk) 

Eote that in equations 4. 11 and 4. 12 the absolute yalue sign has been 


dropped. This is because the summation Indicated will be performed over 


only positive argument and then dbubled. One mus t perform two sums in 
order to determine the value of equation 4. 12 in closed fojrm. They are 
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N-1 

y ~ Z) cos n0 
» = 0 

1^-1 

z = rr S n COS ne 
n«0 

N -1 . 

In 0 

y - re 2- e 
n = 0 


1 . d ^ ine 
^ 'n “ d? ® 

n = 0 


Use of the algebraic identity 1.14 yields 


y = re 


1 - e 
1 - e’ 


d 1 “ e 

z = :;y im 37- — ■ — — 3- 

N d0 . 10 

1 - e 


1 - cos 0 ~ cos W0 + cos (N - 1)0 


2 - cos 0 


1 . r,/ iW0v ,, i0,"2 ^ 10 ie.^ 

= im I (1 - e ) (1 - e ) le - (1 - e ) iNe J 


i 0 \“l _ iWe* 


. / , ilSI0 . TV! ™ 

1 / 1 - e 10 Ne 

- re - " “10 

\ (1 - e ) 1 - e 


. / , fNB -o ^ iN0 

N^®(~ 10.!^ ® (1 " e ) 2 - 2 cos 0 


1 / 1 - e ioV cos N0 - cos (N - 1)0 

N — ^i0^2 ® l ~ 2-2 cos 0 
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1 /I - ie\ . cos (N - 1)0 - cos N9 

z = ® I ^ 2’-"2'^sT'~~ 


11 1-e'^® ie 

y ,z v--^re . g . e 

(1 - e ) 


1 1 


/ 1 


= 2 - N .2+^19 


_ 1 ,1 1 - COS N0 
^ ^ ^ 2 '^N 2-2 COS0 


Substitution of equation 4. 27 into equation 4. 12 gives 


A /U^ + v^ \ I, , 1 /l-cosNA , 1 - COS Ny 
l l I - COST ^ X - ooSY 


(4. 24) 
{4.25) 

(4.26) 

(4.27) 

(4.28) 


Equation 4. 12 for includes terms ranging' from k = 0 to k - N - 1. 

This is only half of the interval desired, since it was already shown that 
terms ranging from -(N - 1) to N - 1 should be included. K one were to 
double a result for equation 4- 12 obtained in equation 4. 28, one would 
make the mistake of including terms with k = 0 two times. As can be 
seen from equation 4. 12, this would give 2 extra contributions, one from 
the term containing A and one from the term containing y* Because 
cos 0 = 1 these two contributions have a total value of 2. Therefore, the 
desired coefficients are 

■■■■■■'■ A' = 2A - 2 (i.29) 

This yields in place of equation 4. 28 
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= 

n 


2 2 / > 
| U V y 1 1 - cos NA ^ 1 “ cos Nv 


2N 


cos A 


1 - cosy 


(4.30) 


It should be noted that this equation is in exact agreement with equation 
3. 16 since ( j c^j ) has twice the power that is actually iri-'each component 
of the frequency spectrum. In addition for equation 3. 16 and 4. 30 to be 
power spectra, they must be divided by odq, the width of a frequency window. 


I V . Symmetry and Spacing 


Since equation 4. 30 was derived for a specific u^ and v , equation 
4. 30 can be written as 


2 2 
u + V 

A- = < ^ . r > 

2N 


n,|j, 


I - cos NA ^ 1 - cos Ny 

I I - cos A 1 - cos Y 


(4.'31) 


or 


A' = 
n 


n,ji 


= S ( 


^ A. 2 

u + V 
M< Y 


1 “ cos N A , 1 - cos N Y 


1 - cos A 


1 - cos Y 


(4. 32) 


This equation is analogous to equation 3 . 19 in the same way that equation 
4.30 is analogous to equation 3.16. 

Equation 4. 31 is a useful starting point for understanding the proper- 
ties resulting from the discrete analysis of a signal which possesses a 
continuous spectrum. This is because of the ease of considering a mono- 
chroir.atic wave, and that results thereby obtained can be linearly super- 
imposed in the limit of a Gaussian time stationary process according to 
equation 4.32. For convenience, rewrite equation 4.31 as 
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1 - cos 2ti (n - u> 

1 - cos 2 ?t {n - p.)/N 


1 - cos 2tt (n + ) 

1 - cos 2 tt (n + li y/K 


where 


e 

u- 



= < 




is proportional to the spec feral power of the g-th component in the continuous 


signal. 


V. Properties of Discretely Sampled Spectra 

Consider the limit in which the \i th component becomes periodic on 
the sampling interval T, that is the limit p, - integer. This is rr -itr • 
sented by 

lim lim ^ 1 - cos 2tt (n - p.) | 1 - cos 2tt (n - u) 

g,-* integer A' =g,-* integer g, 1 - cos 2n(n - g//W g, 1 - cos 2n(n +|j,)/K 

Since real radar systems possess filters which limit the continuous signal 
at the highest unambiguous sample frequency and becatise of arguments 
already made concerning the number of independent Fourier components 
needed for representation of a discretely sampled signal one may write 

: ■ 0 < [j, < r, 0<n<r (4. 35) 

With the restrictions of equation 4. 3 5 
lim 

g, -• integer ^ - 0 


(4. 36) 
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if |i i=n. This is true because the numerators of both terms of equation 
4.34 become 0 in the limit while the denominators do not. If, on the 
.other hand, |i -^ n then 
lim 

|j, -* n AV = (4. 37) 

■■■■■ ■■■■'■ ■■ ' ■ ■' ^ 


The limit is obtained by expanding the numerator and denominator in 
equation 4. 34. Equation 4, 36 and equation 4. 37 can be combined if one 


denotes the limiting process simply by replacing by an integer symbol 


m. This gives 


U + V 

A' =e A = A < ^ 

n, m m n, m n, m 2 . 


Thus equation 4. 34 states that a monochromatic signal which is periodic 

on a sampling interval develops no extra components when sampled. 

This is not surprising since if a signal were periodic on a sample interval 

its nature would be known for all times . Then, in accordance with the 

nncertainiy principle this would mean that the frequency of the mono- 

Ghromatic periodic signal would be completely determined. Another 

statement of the same thing is that periodic signals cannot be spread. 

Now examine equation 4. 34 in the case where tx is not an integer . 

For a particular g , AI has components for all walues of n since 
\h n, |x 

neither numerator nor denominator of equation 4. 34 are zero for any 
Gombination of n and r- . Under these conditions, however; the discretely 
analyzed conbinuous spectrum exhibits certain symmetry properties. It 
has already been stated that in general a particular component is not 
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spread symme trieally aboub Its own frequency. In the special case in 
equation 4.38 the spreading is symmetric since there is no spreading at 
all. There is an interesting non-trivial case, however, in which spreading 
of compoiients not commensui*ate with the sample period is symmetric. 

This case can be understood as follows. In section 2 of the appendix 
it was shown that signals can be analyzed unambiguously from an angular 
frequency of 0 to an angular frequency corresponding to an integer r = 

N/2. The frequency corresponding to the integer r will be called the 
limit of the pass band « Thf center of the pass band is then 

(center of pass band) = ru)Q/2 = TS! u)q/4 (4.39) 

(Note: The numbers in the discrete case corresponding to these values 
are r /2 and N/4. ) 

Consider a case in which there are two components symmetrically 
placed above and below the center of the band pass with equal magnitudes 
and of the form 

s, = u cos 2u (v + N/4)/N + V sin 2tt ('^ + N/4)/N 

and; 

s = u cos 2fr (" V + N/4)/N + v sin 2tt (~ v + N/4)/N (4. 40) 

By equation 4. 33 the components s, and s, possess spectra of the form 
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A' _ 1 - cos 2n (n - v - W4) 

n, VI +N/4 1 - cos 2tt (n - v - TSI/4)/TS! 

1 - cos 2 tt (n + V + N/4) 

- COS 2tt (n + V + N/4)/]SI 

and 

, = 1 - cos 2tt (n - V + lSl/4) 

n, -V + N/4 ® 1 - COS 2it (n - v + N/4)/]S! 


1 " cos 2-rr (n + v - N/4) 

® 1 - cos 2 tt (n + v - ]SI/4)/N 


(4.42) 


To obtain the total spectrum of the signal s, ~ s, + equations 

^tot n -2 

4, 41 and 4. 42 must be added to obtain 


a , ^ r 1 - cos 2 tt (n - V - 1^/4) . 1 - cos 2 tt (n+v ->/4) . 

^tot ^ L 1 - cos 2 tt (n - V - N/4)/]Si 1 - cos 2rr (n -r v - N/4)/]N! ^ 


1 - cos 2tt (n + v t N/4) . 1 - cos 2tt (n - v -t jSl /4) 1 , , 

1 - cos 2tt (n + V + N/4)/lSi cos 2tt (n - v + TSV4 )/nJ 


To see whether a'^ is symmetric about the center of the bandpass 
one must compare with This comparispn 

gives- ' 


Ai - r 1 - cos 2iT(m - V ) 1 - cos 2rr(m ) 

N/4 + LI - cos 2n(m ^ v)/N 1 - cos 2n(m +v )/N 


and 


1 - COS 2Tr(m + v) + Ntt ^ 1 - cos 2rr(m - v ) -f N tt 1 

1 + cbs 2tt (m + V )/N ' 1 + cos 2n (m - v ) /N j 


(4.44) 
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A> = r 1 - cos 2TT(m H-v ) 1 - cos 2TT(m - v) 

N/4 - m, bob ~ ^Ll - cos 2n(m +v )/TSI 1 - cos 2Tr{m - v)/N 


1 - cos 2tt (m - v) + Nit 1 - cos 2rr(m +v) + u 4 

1 + cos 2TT(m - v)/lSl 1 + cos 2n (m + v)/N J 

Since equabion 4. 44 and equabion 4. 45 are idenbical and since superpoSibion 
of ppecbral componenbs holds in bhe bime average iimib ifc has been esbabllshed 
bhat a conbinaous specbrum which is symmebric aboub bhe cenber of bhe pass 
band has bins same symmebiy property when resolved inbo discrebe com- 
ponents. However, in general a conbinuous spectrum is not symmetric 
about the same frequency before and after discrete analysis. The gen- 
eral properfy of asymmetry induced by analyzing a spectrum for a finite 
time is called the principle of spectral asymme try and leads as was 
earlier pointed out to mis-estimates of the spectral center of gravity. 

VI. Spreadihg in the Computer 

Equation 4. 38 states that a signal periodic on the sample interval 
is not spread. Yet in the general discussion of signal and noise cor- 
ruption in the main body of tile report mention was made of corrupting 
only the frequency points which were obtained from a Fourier analysis 
of the backscattered signal. These points correspond exactly to functions 
which are periodic on the sample interval. Taking the autocorrelation 
function of these points triangulating it and re -transforming it to frequency 
space will cause no spreading in the frequency spectrum. One may wonder 
then how, in fact, a frequency spectrum is spread. 
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The spreading is accomplished by interpolation of points between 
tliose actually derived from a Fourier analysis of the incoming signal. 

If N samples are taken on a sample interval ISI autocorrelation functions 
and N/2 frequency points result. Before re transforming to autocorrela- 
tion space, ISI / 2 frequency points are added. These points are interpola- 
tions between the N/2"real" points already obtained. Wow there are N 
points in frequency space. These W points are not all independent, since 
W/ 2 points are simply interpolated from the "real" points. The N points 
are symmetrized and the Fourier cosine transform is taken. That is, 

2W points of a symmetric interpolated spectrum are transformed into a 
2N point autocorrelation function. This function is symmetric and so W 
points are unnecessary since they contain no new information. Thus, 
only N points are used and these points are the new W point autocorrela- 
tion function, which is triangulated to produce spreading in frequency 
space. The interpolated points in frequency space are not periodic on 
the sample interval and thus triangulation causes these points to be spread 

To examine whether or not the spreading is of the proper magnitude, 
that is whether it has the same value as it would in a continuous spectrum, 
the following argument is employed. Imagine a continuous spectrum ex- 
tending between two "real" (uninterpolated) points in frequency space. 
Assume that the interval under consideration contains one interpolated 
and one uninterpolated point. This interval Is therefore representative 
of the way in which points were added in frequency space . For conven- 
ience allow the continuous spec trum to have an amplitude 1. Assume 
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the spectrum has height 0 everywhere except between tlie m and (m + l)th 
real frequency points, with the (m -I- l)th frequency point being excluded. 
By equation 4. 32 the discretely analyzed spectrum has the form 


- cos 2n(n - |x) ^ 1 - cos (n +ii.) 

cos 2Tr(n - 1 ■ cos 2Tr(n +|i)/IS), 


The denominators are slowly changing functions by comparison to the 
numerators of equation 4.46. Thus, over the small interval m, m + 1 
to a fii’st approximation the denominators of the function of equation 4. 46 
can be considered constant and evaluated at m + I while the numerators 
can be integrated. This yields 


A' J.1 
n, m + 2 


cos 2n(n - [m + 1/2] )/N i - cos 2 it (n + M + 1 / 2 ) /N 


{4.47) 


[Since the cosine terms in the numerator cycle through one period between 
m and m + 1, hence integrate to 0. Now consider the discrete analog of 
equation 4. 47. First in equation 4. 47 the spectrum wHch was used had 
an amplitude pf 1 over an interval of length 1. Thus the total pove r in 
the frequency interval between m and m + l was unity. To make a con- 
sistent analog in discrete space the total power must be the same. Thus, 
the powers in the points at m and m + 1 / 2 must be equal with amplitudes 
1 /2. It has already been shown that the point at m will not be spread. 
Thus, the spread power in the spectrum comes from the point at m + 1/2, 
and lay equation 4. 32 
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A . _ m + 1/2 _ 1 / 2 

n “ 2 ■ 2 VI - cos 2n [n - (m + 1/2)] /N 

■■■ "'v ' 2 ' "" ' ■ 

1 - cos 2n [n + (m + l/2)/N] 

Equations 4.47 and 4.48 are identical. Of course, they are only 
approximately so since equation 4. 47 was obtained assuming that the 
denominators in equation 4. 46 were constant. However, this is a good 
approximatibn except where these denominators are near 0. In such a 
case, however, components nearly pei’iodic on the sample interval are 
being considered. As the denominator passes through 0 there is a positive 
and a negative contribution to the integral and these contributions tend 
to cancel to 0. This cancellation simply confirms the fact that nearly 
periodic functions tend to have little spreading. (Note: The same results 
for spreading are obtained if more than one point is interpolated for each 
real point. In fact, the approximate result obtained in equation 4. 47 
becomes exact in the infinite limit. ) 


(4.48) 


